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B s — > /o(980) form factors and B s decays into / (980) 



O 

(N 

PL. 

in 



Oh' 
i ' 

D 
43 



> 
O 
oo 

00 
(N 

(N 

o 
o 



X 



Pietro Colangelo, Fulvia De Fazio and Wei Wang 
Istituto Nazionale di Fisica Nucleare, Sezione di Bari, Italy 

We compute the B s — > /o(980) transition form factors using light-cone QCD sum rules at lead- 
ing order in the strong coupling constant, and also including an estimate of next-to-leading or- 
der corrections. We use the results to predict the branching fractions of the rare decay modes 
B s ->■ f £ + £~ and B s -> f uu, which turn out to be O(10" 7 ) (B s ->■ f (980)£ + r, with I = e,fi), 
0(1CT 8 ) (B s -> /o(980)t + t~) and 0(1CT 6 )(B 3 -» /o(980)i/p). We also predict the branching ratio 
of B s — s> J/^/o(980) decay under the factorization assumption, and discuss the role of this channel 
for the determination of the B 3 mixing phase compared to the golden mode B s —¥ J/ip(f>. As a last 
application, we consider D s —¥ fo form factors, providing a determination of the branching ratio of 
D s -¥ f e + u e . 

PACS numbers: 13.20.He, 13.20.Fc, 13.25.Hw 



I. INTRODUCTION 



Theoretical and experimental efforts aimed at disclos- 
ing physics beyond the standard model (SM) proceed in 
several directions. Among these, there is the study of 
rare processes which are induced only at loop level in 
the SM and are therefore sensitive to new physics (NP) 
contributions which may potentially enhance their small 
(< 10~ 5 ) branching ratios Another testing ground is 
the precise study of CP violation. It has been realized 
that the amount of CP violation within the SM is too 
small to explain the observed baryon asymmetry of the 
Universe [2|, a conclusion confirmed by recent analyses 
||. Since the only source of CP violation in the SM is 
the complex phase of the Cabibbo Kobayashi Maskawa 
(CKM) mixing matrix, the determination of the elements 
of this matrix and of their relative phases is of primary 
importance, in order to disentangle souces of additional 
contributions to CP violation. As well known, the task 
is afforded through the study of the so called unitarity 
triangles, the graphical representations of the conditions 
stemming from unitarity of the CKM matrix. The most 
studied triangle is the one which relates the CKM ele- 
ments involved in B decays. Direct and indirect deter- 
minations of its sides and angles lead to a picture of CP 
violation coherent with the SM description. Also in this 
case, investigation of effects predicted to be small in the 
SM is a promising strategy to reveal new physics. 

B s mesons provide the possibility to search for new 
physics scenarios exploiting both the strategies outlined 
above. On the one hand, rare B s decays induced by b — > s 
transition are suppressed in the SM, as with all decay 
modes governed by such a transition, and new physics 
effects may enhance their branching fractions. For ex- 
ample, it has been shown that, in presence of a single 
universal extra dimension compactified on a circle with 
radius R, the rates of B s — > (f>vv, B s -> ifk+£- and 



B K 



are enhanced when R decreases, while 



the opposite happens in the case of B s — > (fry, which has 
a smaller branching fraction with respect to SM for small 
values of R" 1 [|. 



On the other hand, the analysis of the unitarity tri- 
angle of CKM elements relevant for B s decays is an im- 
portant test of the SM description of CP violation. The 
triangle is defined by the relation 



v us v: b + v cs v; b + v ts v t t = . 



(i) 



Vc,V*, 



is 



One of its angles, j3 s , defined as f3 s = Arg 

half of the phase of B s — B s mixing, and is predicted to 
be tiny in the SM: /3 S ~ 0.019 rad. Recent data obtained 
by the CDF @ and DO Q Collaborations, based on the 
angular analysis of B s — > J/ip<p, indicate much larger 
values, although with sizable uncertainties, so that the 
precise measurement of /3 S represents one of the priorities 
in the physics programs at the hadron colliders and at the 
B factories operating at the T(5S) peak 0]. 

In this paper we consider B s decays in both respects. 
We compute the B s — > /o(980) 1 form factors using light- 
cone QCD sum rules (LCSR) at the leading order in the 
strong coupling constant (Sect. IIIIIII Aj) and including 
an estimate of next-to-leading (NLO) corrections (Sect. 
IIIIB[) . In Sect. IIV Al we use the results to predict the 
branching fractions of the rare decay modes B s — > fo£ + £~ 
and B s —> fovv in the SM. The form factors are also 
a necessary ingredient to study the nonlcptonic mode 
B s —t J/tpfo which, together with B s —> J/tp(j), permits 
one to access the phase f3 s Q. Our predictions for this 
mode are collected in Sect. IIVBI As a byproduct of the 
calculation, we explore the D s — > f$e + v e decay channel, 
the branching ratio of which has been recently measured 
by the CLEO Collaboration @,[l(j. Conclusions are pre- 
sented in the last section. 



1 Hereafter, we use fo to denote the /o(980) meson. 
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II. LIGHT-CONE QCD SUM RULE 
CALCULATION OF B s -> fo FORM FACTORS 

The matrix elements involved in B s — > fo transitions 
can be parameterized in terms of form factors as 



The current jr 2 interpolates the B s meson: we choose 
jr 2 = bij^s. Its matrix element between the vacuum and 
B s is given in terms of the decay constant fs a , 



(B s ( PBs )\bi l5 s\0) 



-f Bs - 



(7) 



(fo(Pf )\s-f^ 5 b\B s { PBs )) = 



Tor — mi 

D. JO 



+ F (q 2 )- 



J i> 



(fo (Pfo )\s<?ni>ls<fb\B s (p Bs )) 
F T (q 2 ) 



"IB, 



>J\, 



)q,, 



(3) 



where P = p B , + Pf and q = p Bs — Pf - In this section 
we describe the calculation of the three functions Fx , F 
and Ft using the method of light-cone QCD sum rules. 
For the sake of the calculation, it is convenient to define 
the auxiliary form factors /+ and /_, 

(h(Ph)\sWb\B s (p Ba )) = -l{.f+{q 2 )P» + f-{q 2 )q»} 

(4) 

in terms of which Fx and Fq read: 



Fl(q 2 ) = f+(q 2 ) , 
Fo(q 2 ) = f+(q 2 ) + 



2 2 

m „ — mi 

i>s Jo 



f-(q 2 ) 



(5) 



As a reconciliation of the original QCD sum rule ap- 
proach [ll[ and the application of perturbation theory 
to hard processes, LCSR [l2j present several advantages 
in the calculation of quantities such as the heavy-to-light 
meson form factors. The method includes both hard 
scattering and soft contributions. In the hard scatter- 
ing region the operator product expansion (OPE) near 
the light-cone is applicable. Based on the light-cone 
OPE, hadronic quantities, like form factors, are expressed 
as a convolution of light-cone distribution amplitudes 
(LCDA) with a perturbatively calculable hard kernel. 
The leading twist and a few subleading twist LCDA give 
the dominant contribution, while higher twist terms are 
power suppressed. The LCSR approach has been suc- 
cessfully applied to compute the hadronic parameters in- 
volved in many different processes fl3j . 

The starting point for a LCRS evaluation of form fac- 
tors is the correlation function of suitably chosen quark 
currents. Here we consider the correlation function 

Tl(p f0 ,q) = i [ d 4 xe^(fo(p fo )\T{ Jri (x),j r2 (0)}\0) 



(0) 

where jr 1 is one of the currents appearing in the matrix 
elements (J2J3J) defining the form factors: jr 1 = Jfj = 
s7a<75& for F x and F 0l and j Tl = J^ T = scr^j 5 q L 'b for F T . 



(2) 



TOfe + TO, 

We also introduce the /o(980) decay constant // , 

(/oO/ o )M0) = to/o//o (8) 



needed in the following; ff has been evaluated by several 
groups [T3 - [T7i . 

The LCSR method consists in evaluating the corre- 
lation function Eq.© both at the hadronic level and 
in QCD. Equating the two representations provides one 
with a sum rule suitable to derive the form factors. 

The hadronic representation of the correlation function 
in © 

n HAD (p f0 ,q) = 

(h{Ph)\jrAB s {p h +g))(B s (p fo +g)|jr 2 |0) 

m B 3 - (pfo + q) 2 

(/o(P/o)b"ri \HPfo + q))(h(Pf + g)|jr 2 |0) 



E 



m l - (pa + q) 2 



(9) 



consists in the contribution of the B s meson and of the 
higher resonances and the continuum of states h. In 
a one-resonance+continuum representation, the correla- 
tion function can be written as 



n 



HAD 



(Ph > q) 



(h{Pfo)\kAB s {p h +q))(B s (2> fa +q)\jrM 



- (Ph + qf 



ds- 



P h (^q 2 ) 
- (Pfo + q) 2 ' 



(10) 



where higher resonances and the continuum of states 
are described in terms of the spectral function p h (s,q 2 ), 
which contributes starting from a threshold sq- 

At the quark level, the correlation function can be eval- 
uated in QCD in the deep Euclidean region, where it can 
be written as 



1 [°° ImnQ CD (.s,g 2 ) 

-/ ds —^.(11) 

77 J( mb+ms y- s -(Pf +q) 



n QCD (p /o ,g) 



Using global quark-hadron duality, the integral in (|10p 
can be identified with the corresponding quantity in the 
QCD representation (fTTj) : 



ds- 



P h (s,q 2 ) = i 
- (Ph + q) 2 k 



ds 



lxxxU^ CB (s,g 2 ) 

s - (Ph + q) 2 



(12) 



A Borcl transformation of the hadronic and of the QCD 
expressions of the correlation function is carried out, de- 
fined as: 



S[J(Q 2 )] = lira 
1 



2\n 



(n-1)! 



(-Q 2 ) 



oo, n— too. 

d 



dQ 2 



9±=M* 

HQ 2 ) , 



(13) 



where 3" is a function of Q 2 = —q 2 and M 2 is the Borel 
parameter, so that 



B 



1 



(s + Q 2 )" 



exp(-s/M 2 ) 
(Af 2 )« (n-1)! 



(14) 



This operation improves the convergence of the OPE se- 
ries by factorials of n and, for suitably chosen values of 
M 2 , enhances the contribution of the low lying states to 
the correlation function. 

Applying the transformation to both n HAD and n^ CD 
we obtain the sum rule 



,2 n 



(/o (p/o )\3Ti\B. (pu. )) (B s ( PBs ) \jr 2 1 0)cx P 

dsexp[-s/M 2 ] Imn QCD (s,q 2 ), 



M 2 



(15) 



" ^ (ms+m, 

where ps a = Pf + q- Ea.(|15| allows to derive the sum 



rules for /+, /_ and Ft, choosing either the current jr 1 = 
or the current jn = J^ T ■ 

The calculation of n^ CD is based on the expansion of 
the T-product in ^ near the light-cone, which produces 
matrix elements of non-local quark-gluon operators. In 
the description of fin as a ss state modified by some 
hadronic dressing [14} . these can be defined in terms of 
fa light-cone distribution amplitudes of increasing twist: 



</ (p/oMa07^(0)|0) = ff Pf » / due™***-** fQ (u 



(/o(p/ )|*(xMO)|0) = m /o / /o / dtie^/o-^^), 

Jo 

Tfi f — 

(/o(p/ )|s(a;)^ s ( )l ) = — ^//oIP/o^-J'/o^m) 

x / dtte <, *A> (16) 
Jo 

where the LCD A $/ is twist-2, and the other two are 
twist-3, and are normalized as 

du# /o (u) = 0, f du$%(u) = f du<f>%(u) = 1. 
io Jo Jo 

(17) 

In terms of these LCDA, the sum rules for the three 
form factors read: 



J 



F T (q 2 ) 



nib + m s 
2m 2 f Bs 



// exp 



' m B s ' 


f r 1 du 


m 2 + uurn* — uq 2 


( / — exp 

[J u U 


. M 2 _ 


uM 2 



-m b $ fo (u) + um fo <5> s h {u) + -m h § a fo (u) 



ml + q 2 



u 2 m 2 fo m fo $ a fo (u) 



uM 2 



■ exp [s /M 2 } 



// exp 



~ m B s ~ 


\ f 1 du 


m 2 + uum 2 o — uq 2 


< / —exp 

[J u B U 


. M2 . 


uM 2 



171b + TUg 

2m 2 B J Bs 

1 — u „ , 

uo{m 2 + q 2 - ulmjj + 2 (to 2 - q 2 + u 2 m 2 fo ) 
u {m 2 + ulm 2 fo -q 2 ) 

(mB s + mf )(mb + m s ) 



6 to 2 + MgTO 2 o — q 2 



[m b $ h (u) + (2-u)m fo $%(u) 



(18) 



u(m 2 + q — u 2 m 2 o ) + 2(m 2 — q- + u 2 m^ o ) to/ $j q (w) 



u 2 M 2 

2^,2 



"//o ex P 





f fi 1 dw 


< / —exp 


. M2 . 



r /,r2 1 TO /o^/o("o) 

exp [—s /M\ 



(to 2 — uq 2 + uum 2 Q ) 
uM 2 



(19) 



h mb- 



6uAf 2 



-m b 



"Vo$/ M exp[-s /M 2 } 
6 to 2 — q 2 + u^m 2 a 



(20) 



where 
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Our formulae can be compared to the ones for the B-to 
scalar meson form factors in Ref. 
the meson ao is considered. We find differences in the 
expression of the form factor / + . 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

A. Leading order results 

Based on the conformal spin invariance, the LCDA can 

3/2 

be expanded in terms of Gcgenbaucr polynomials C n ■ 
The expansion of the twist-2 LCDA <J>/ (u) reads: 

<P fo (u)=6u(l-u)^B + J2B n Cf/ 2 (2u-l)^ . (22) 

Due to the charge conjugation invariance, all even Gegen- 
bauer moments of $ f (u) vanish, so that B2 m = for 
to = 0, 1, ■ ■ ■ in (|22]k as for the odd moments, we in- 
clude only the first one, using the value of the coefficient 
Bi = -0.78±0.08 fixed in ref. 0. For the twist-3 LCDA, 
due to the lack of knowledge about their moments, we use 
the asymptotic form, i.e. the first term of the Gegenbauer 
expansion, 



results for Fi(0), a, and bi are collected in Table U and 
18j , where the case of the q 2 dependence is depicted in Fig. [2] The uncertain- 



= 1, *5 b (u)=6u(l-«). 



(23) 



Let us quote the numerical values of the other physical 
parameters. The meson masses are fixed to the PDG val- 
ues rriB s = 5.366 GcV and m/ = 0.98 GeV [lj|, while for 
quark masses we use nib = 4.8 GcV and m s = 0.14 GeV 
[m I20I ]. As for the decay constants, we use Jb, = 
(0.231 ± 0.015) GeV [2l| and f fo = (0.18 ± 0.015) GeV 
[Hj] 2 . The threshold s is fixed to s = (34 ± 2) GcV 2 , 
which should correspond to the mass squared of the first 
radial excitation of B s . 

With these numerical inputs, the sum rules (fl"8T) -([20" |) 
provide us with the form factors for each value of q 2 as a 
function of the Borel parameter. The result is obtained 
requiring stability against variations of M 2 . 

In Fig. [1] we show the dependence of the form factors at 
q 2 = on the Borcl parameter M 2 . We observe stability 
when M 2 > 6 GcV 2 , and we fix M 2 = (8 ± 2) GeV 2 . 
To describe the form factors in the whole kinematically 
accessible q 2 region, we adopt the parameterization 



F.iq 2 ) = 



Fi(0) 



1 - aiq 2 /m% + bi{q 2 /m 2 B 



(24) 



where Fi denotes any function among -F^o.t- The param- 
eters di , bi are obtained through fitting the form factors 
in the small q 2 region (we choose < q 2 < 15 GeV 2 ); the 



ties in the results reflect those of the input parameters. 
In Table U we also report the values of the form factors at 
zero-recoil (g 2 lax ) which are derived using the expression 
in Eq. (|23)l. 

The results in Tabic U show that the parameters a; and 
bi determining the q 2 dependence are close to each other 
in the case of F± and Ft- The reason is the following. In 
the heavy-quark limit and in the large energy (LE) limit 
of the recoiled meson, the three B s — > fo form factors 
can be related to a single universal function £/ which is 
specific for fo and does not depend on the Dirac structure 
of the current appearing in the definition of the various 
matrix elements, such as those in Eqs. (J2H]) [22}. When 



the energy E of the light meson in the the final state is 
large, such relations read as 



m Bs + to/ 



F T (q 2 ) = F^q 2 ) 



2E 



W), (25) 



where, neglecting to 2 o but keeping to/ in the kinematical 
factors, E is related to q 2 : 



- 2m Bs E. 



(26) 



The first equality in Eq. ([25]) shows that the large en- 
ergy limit predicts that F\ and Ft have the same q 2 
dependence. For the shape parameters of Fq, one can 
obtain two relations through the second equality: 



a = — l + oi, 60 = 1 — 01+61. 



(27) 



0.6 
0.5 
0.4 

F,(0) 0.3 
0.2 
0.1 

0.(1 



0.5 
0.4 

F r (0) 0.3 
0.2 
0.1 

0.(1 



M 2 (GeV 2 ) 



In ref. Qj a larger result is reported: f fo = (0.37 ± 0.02) GeV. 



M 2 (GeV 2 ) 



Figure 1: Dependence on the Borel parameter M 2 of the 
B s — > fo form factors at q 2 = 0: -Fi(O) = _Fo(0) (upper panel) 
and Ft(0) (lower panel). 
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Table I: Parameters of the B s —¥ fo form factors by LCSR 
at the leading order. The values of F i (cfi n ^) are evaluated 
through Eq.fUJl. 



Fi(g 2 = 0) en 



In 



Fi{q 



'max; 



Fi 0.185 ±0.029 1.44±g;J| 0.59 
Fo 0.185 ±0.029 OAlt^ol O.OltHg 
F T 0.228 ±0.036 1.42±S ^ 0.60 



,+0-07 n m4 +0.158 
0.05 u - ul *-0.102 
+0.055 
0.038 

,+0.06 n 71 ..+0.197 
-0.05 u -' lq: -0.126 



Table II: B a — >• /o(980) transition form factors obtained in- 
cluding an estimate of next-to-leading order corrections (see 
text). 



Fi(<7 2 = 0) 



b, 



+0.13 



Fi 0.238 ±0.036 1.50 
F 0.238 ± 0.036 0.53i[J;io 
Ft 0.308 ±0.049 lA6±° a ]i 



0.58 



+0.09 
-0.07 

-0 Sfi+ 09 
u-JO-o.os 



0.58 



+0.09 
0.07 



1.0 rr 
0.8 - 



0.6 - 
0.4 - 
0.2 - 




0.0 L — ' — ■ — ' — ' — 1 — 1 — 1 — ■ — ' — 1 — ■ — 1 — ■ — ■ — L - 
5 10 15 

q 2 (GeV 2 ) 

1.0 p , , , , 1 , , , , 1 , , , , r- 



0.8 



0.6 

Foiq 2 ) 

0.4 



0.2 




0.0 L — ' — ■ — ' — ' — 1 — ' — 1 — ■ — ' — 1 — ■ — 1 — ■ — ■ — L - 
5 10 15 

^(GeV 2 ) 

1.0 p . . ■ ■ 1 . ■ ■ . 1 ■ . . . r- 



0.8 - 
0.6 - 

F T (q 2 ) 

0.4 - 

0.2 - 




0.0 b — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — J 
5 10 15 

g 2 (GeV 2 ) 

Figure 2: q 2 dependence of the B a — > fo form factors. 

Using the results for a\ and b\, we find from (f27|): 

a[ LE) ~ 0.44±0.1 and b[ LE) ~ 0.15 ±0.12; therefore, the 
first relation in (|27| is well respected in our calculation, 
while not much can be said about the second relation due 
to the uncertainty affecting bo- 



B. Estimate of the next-to-leading order 
corrections 

In order to provide an estimate of next-to-leading or- 
der effects in the determination of the B s — > fo form fac- 
tors, it is worth comparing this case to the calculation of 
B — > 7r form factors. In B — > ir transition, both the light 
quarks and the light 7r meson have small masses which 
can be safely neglected, while the strange quark and the 
scalar meson fo masses may induce sizable effects. An- 
other observation is that, neglecting the quark masses, 
the Lorentz structures of pion and fo matrix elements 
differ by a minus sign in terms proportional to the twist- 
2 LCDA. Finally, contributions from the twist-3 LCDAs 
in B — > 7T transition are characterized by the chiral scale 
parameter fi v , while in B s — > fo they arc proportional to 
the mass of fo- 

In LCSR, NLO corrections to B — > tt form factors have 
been studied by two groups [H, H3] , while the complete 
expressions for the NLO corrections to B s — > fo form fac- 
tors are not known at present. The expressions relevant 
for B — >• 7T form factors given in Ref. |24j can be used to 
estimate the radiative corrections in the case B s — > fo, 
keeping in mind the three differences above. We first con- 
sider the changes to the leading order result due to the 
different treatment of quark and hadron masses. Setting 
the quark mass m s to zero, the values of the form factors 
are reduced by about 3%. The mass of fo, the analogous 
of the pion mass m n and the chiral scale parameter 
cannot be put to zero, as this would smear all terms from 
twist-3 LCDA: we set the mass square of fo to be zero 
keeping the linear terms in the form factors, obtaining 
an enhancement of the form factors by about 3%. Af- 
ter that, evolving all the scale-dependent parameters to 
a scale of about the Borel mass, /j, ~ 3 GeV, we find that 
the leading order contributions are further ly enhanced, 
obtaining the central values: Fi(0) = F o (0) = 0.216, 
oi = 1.50, h = 0.58, a = 0.216, b = 0.53 and 
F T (0) = 0.262, a t = 1.46, b t = 0.58. Then, radiative cor- 
rections to twist-2 and twist-3 LCDA are also found to 
be rather small, the B s — > fo form factors being changed 
to Fi(0) = F (0) = 0.238 and F T (0) = 0.308. The result- 
ing values, with the inclusion the uncertainty due to the 
input parameters, are collected in Table [TTJ they are also 
used in the phenomenological analysis, keeping in mind, 
however, that the procedure used in their determination 
must be considered as only approximate. 

Before closing this section, it is worth mentioning 
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Table III: B s ->• /o(980) form factors at q 2 = 0. Results 
evaluated by CLFD/DR [H, PQCD [H and QCDSR [13] 
approaches are collected for a comparison. 



quark and gluon fields: 



CLFD/DR 


PQCD 


QCDSR 


This work 


Fi(0) 0.40/0.29 a 
F T (0) 


n oc+0.09 b 
n 4 n+0.10 6 


0.12 ±0.03 c 
-0.08 ± 0.02 c 


0.185 ±0.029 
0.228 ±0.036 



"using f_ Bs = 0.259 GeV 

6 using f Jo = 0.37 GeV 

c using f fo = 0.37 GeV and f Bs 



0.209 GeV. 



that the B s — > /o(980) form factors have been com- 
puted by other approaches: the method based on covari- 
ant li ght- front dynamics (CLFD) and dispersion relation 
(DR) [HI, the perturbative QCD approach (PQCD) [H, 
short-distance QCD sum rules (QCDSR) [27|. The re- 
sults are collected in Table IIIII The form factors by 
PQCD are proportional to the /o decay constant, while 
those by short-distance QCD sum rules are proportional 
to the inverse of this constant. Thus, a larger decay con- 
stant, // = 0.37 GeV as reported and used in gives 
larger form factors in the PQCD approach and smaller 
ones in QCDSR with respect to ours. Taking into ac- 
count the difference in the decay constant, the results 
in Refs. (2(1 I27J are consistent with ours, while the two 
results in Ref. (25| are sensibly larger. 



IV. PHENOMENOLOGICAL APPLICATIONS 



A. Semileptonic B s — > fo£ + £ and B s — > favv decays 



As a first application of our study, we predict the 
branching ratios of the decays B s — > fo£ + £~ and B s — > 
Jqvv, processes which, being induced by the flavor- 
changing neutral current transition b — ¥ s, are potentially 
important for detecting new physics effects. 

The SM AB = 1, AS = — 1 effective Hamiltonian 
describing the transition b — > s£ + £~ can be expressed in 
terms of a set of local operators: 



r 10 



(28) 



Gf = 1.166 x 10 5 GeV 2 being the Fermi constant and 
Vij the elements of the CKM mixing matrix (since the 

is O(10~ 2 ), we neglect terms proportional 

The operators Oi are written in terms of 



rati0 V V* 

Vtb Vts 

to V uh V* s ) 



01 = 

2 = 

3 = 

Oi = 

o 5 = 

o e = 

o 7 = 

Os = 

9 = 

010 = 



(sLa7 A1& L^)(cL/37M C ia) , 
(sLa^b La )[(u Ll 3'Y tl U L p) 



{s La l^b L p)[{uLpl^u La ) 

(sLc7 A1& ia)[(wi?^7M U -R/3) 
{mbSLaO-^bRa 

~2 



16tt 2 

16tt 2 

e 2 
16tt 2 

c 2 



• ± {bLDl^b L f})] , 

• ± (pLPl^Lot)] , 

■ + (bRplnbiip)] , 

■ + (bRplnblla)] , 

m s s Ra <r^b La ) , 



SLaV 



j.i v 



^r) b 



R;3 



(29) 



1±75 



6, and cr^ 



16tt 2 

with a, f3 color indices, bn,L 
i 

— [7^, 7"]; e and g s are the electromagnetic and the strong 

coupling constant, respectively, and and G^ v in O7 
and Os denote the electromagnetic and the gluonic field 
strength tensor. 0\ and O2 are current-current opera- 
tors, O3, Oq QCD penguin operators, O7 and Os mag- 
netic penguin operators, Og and O w semileptonic elec- 
troweak penguin operators. The Wilson coefficients in 
are known at NNLO in the Standard Model (281. 



The operators Ox and O2 contribute to the the final state 
with a lepton pair through a cc contribution that can give 
rise to charmonium resonances J/ip, ip(2S), • • • , resonant 
term which can be subtracted by appropriate kincmati- 
cal cuts around the resonance masses. The Wilson coef- 
ficients C3 — Cg are small, hence the contribution of only 
the operators O7, Og and Oio can be kept for the descrip- 
tion of the b — > s£ + £~~ transition. In our study we use 
a modification of the Wilson coefficient C7: , which 
is a renormalization scheme independent combination of 
Ct,Cs and C2, given by a formula that can be found, 
e.g., in P. 

The B s and fa matrix elements of the operators in 
([29]) can be written in terms of form factors, so that the 
differential decay width of B s — > fo£ + £~ reads: 

dT(B s -> fot+t-) 



dq 2 



G 2 F a 2 em \V tb \ 2 \V t * s \ 2 V\ q 2 -4m 2 1 



512m| s 7r 5 



V q 2 3q 2 
6m 2 \C w \ 2 (m 2 Bs ~m 2 fo ) 2 F 2 (q 2 ) 



(30) 



±(g 2 ± 2m 2 e )X 



2 , 2C 7 {m b - m s )F T {q 2 ) 



rriB s + mf 



+ \C w \ 2 {q 2 -Ami)\Ff{q 2 ) 
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JO' 



■4m 2 fo q 2 , 



with A = A(m| a , m\ , q 2 ) = (m% s - q' 
a em = 1/137 the fine structure constant and mi the lep- 
ton mass. 

Analogously, the SM effective Hamiltonian for b 
svv. 



V tb V* Vx X(x t )0 L = C L O l 



(31) 



a/2 2tt sin 2 ( 
includes the operator 

O l = - 75 )6) (P7 M (1 - 76)1/) . (32) 

o 

771 

6*vf is the Weinberg angle; the function X(xt) (xt = — j-, 



with m t the top quark mass and the W mass) has 
been computed in [3(1 and (UH^I, while the QCD factor 
rjx is close to one [3ll433l |. so that one can use rjx = 1- 
From this effective Hamiltonian, the differential decay 
width 



dT(B s -j. f uD) 
dq 2 



|C L | 2 A 3 / 2 (m|,,m 2 /o ,g 2 ) 



2 \l 2 



96m 3 3 7r 3 



can be obtained. 

In the numerical calculation we use 

C 7 = -0.30137, C* 9 = 4.1696, C w = -4.46418, 

C L = 2.62 x 10~ 9 , (33) 

together with V ts = 0.0387 and Va = 0.999 0. Using 
these inputs and t(B s ) = 1.47 ps [19| we find: 



B3l(B, -> faTtT) = (9.51 



+3.1^ 



x 10" 



Sft(5 s -> /or+r") = (l.l±g:|) x 10 



(34) 



Sft(5 s / i/i7) 



-T±|;|) x 10' 



with £ = e,fi. Our estimate of the NLO effects in the 
form factors modifies the branching ratios to CBfR(i3 s — > 
/ £+£-) = (16.7 ±6.1) x 10" 8 , Bft(i? s -> /or+r") = 
(2.7 ± 1.3) x 10~ 8 , and S3?(B S -> / z/i/) = (15.2 ± 5.6) x 
10~ 7 . These decay modes arc therefore accessible at the 
LHCb experiment at the CERN Large Hadron Collider 
and at a Super B factory operating at the T(5S) peak. 



B. Nonleptonic B a — >■ J/xpfo transition 

The study of CP violation and the measurement of 
the CKM angles mainly proceed through the measure- 
ment of nonleptonic decay modes. In the B s sector, 
the channel B s — > J/ip<p is the golden mode to inves- 
tigate CP violation, and from the analysis of this mode 
the CDF 01 and DO @ Collaborations at the Fermilab 
Tevatron have obtained values of the B s mixing phase 
4> s = —2(3 S much larger than predicted in the SM, mod- 
ulo a large experimental uncertainty. If confirmed, this 
measurement would indicate physics beyond SM. It is of 



prime importance to consider other processes allowing to 
access /3 S , namely B s — > J/tpf], Jji\)i]' and J/i/>/ (980) 
in which the final state is a CP eigenstate and no an- 
gular analysis is required to disentangle the various CP 
components, as needed for B s — > J/ip(j). However, the re- 
construction of B s modes into 77 and 77' is experimentally 
challenging, since the subsequent rj or 77' decays involve 
photons in the final state. The case of /o seems feasible, 
since /o(980) essentially decays to tt + tt~ and to 2tt° (the 
decay to K + K~ has also been seen) [lj|. Theoretical 
predictions of B s — > J/tpj b(980) are therefore of great 
importance. 

The quantitative description of nonleptonic decays is 
very challenging. The theoretical framework to study 
such decays is based on the operator product expansion 
and renormalization group methods, which allow to write 
an effective hamiltonian as in the case of the modes con- 
sidered in the previous section. However, now one has to 
consider hadronic matrix elements (J/ipfo\Oi\B s ) with 
Oi four-quark operators, the calculation of which is a 
nontrivial task. One of the strategies which has been ex- 
ploited is the naive factorization [34j | , in which such quan- 
tities are replaced by products of matrix elements of the 
weak currents appearing in each one of the operators of 
the effective hamiltonian relative to the considered pro- 
cess. These objects are expressed in terms of meson decay 
constants and hadronic form factors. This procedure is 
affected by several drawbacks, and various refinements 
have been proposed. It has been shown that a theoret- 
ical justification of naive factorization in the case of B 
decays can be found in the heavy quark limit mj, — > 00 
only in a limited class of processes [35j. One can con- 



sider the so called generalized factorization approach, in 
which the Wilson coefficients (or appropriate combina- 
tions of them) appearing in the factorized amplitudes are 
regarded as effective parameters to be fixed from experi- 
ment, a procedure adopted in the following. 

Using the factorization ansatz, the decay amplitude 
B s {pb b ) ~> J/tp(p,l>,e)fo(Pf ) (e being the J/ip polar- 
ization vector, pb s , Pin Pf the momenta of the three 
particles) is given as 

G F 

A{B S -> J/ipfo) = -^=VcbV* s a%m^fji^ 



(m 2 /^* -psj; (35) 



/ju is the J/ip decay constant, determined from the 
J/tp ->• e+e' decay width 0: / J/v , = (416.3 ± 5.3) 
MeV. The factor 02 is a combination of Wilson coeffi- 
cients which can be extracted from the B — > J/ipK de- 
cays, under the assumption that is the same in the 
two processes. For these decays the branching ratios are 
known 11911: 



•BJl(B- - 
T,3l(B 



J/ipK~) = (1.007 ± 0.035) x 10" 3 , 
>J/ipK°) = (8.71 ± 0.32) x 10~ 4 . (36) 

In order to extract a-i-, the form factor F^~^ K is re- 
quired. We use two different paramctcrizations, obtained 
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by short-distance (CDSS) [36| and light-cone QCD sum 
rules (BZ) [23}. The result is different for the two sets of 
form factors, while there is almost no difference whether 
we use the charged or the neutral channel: 

iB->J/iJ>K,(CDSS) I n o QzL +0.053 
\ a 2 I — u - oaq: -0.04lJ 

| a s^j/^(sz)| = 25±0 03 (37) 

To be conservative with the hadronic uncertainty, we use 
the average value 0,2 = 0.32 ± 0.11 of the two values to 
compute T>3l(B s J/^/o)- Using V cb = 0.0412, V cs = 
0.997 [lj| and our LO prediction for the B s — > f form 
factors, we obtain 

^{B s J/^/n) = (3.1 ± 2.4) x 10" 4 (38) 



while, including our estimate of NLO corrections, the 
branching fraction is T>Jl(B s J/tpfo) = (5.3 ± 3.9) x 
10~ 4 . The rate of B s —> J/ipfo is large enough to per- 
mit a measurement; notice that the branching fraction of 
B s -> J/ip(j) is -BJi(B s -> J/ip4>) = (1.3±0.4) x 10- 3 [H. 



To gain a better insight on this point, it is interest- 
ing to compare these results to the branching fraction of 
B s — > J/t/jl4>l {L denotes a longitudinally polarized me- 
son) computed in the factorization approach. Neglecting 
the mass difference between <p and fo in the phase space, 
the ratio of the branching fractions of the two modes can 
be written in terms of form factor combinations: 



R 



■BR(B S -> J/V/o) 



[Ft 



(ml)] 2 \(m 2 Bs ,ml,m 2 fo ) 



[A? 



(m? 



2m d 



- Ai 



Mi 



0.13 ±0.06 
0.22 ±0.10 



(39) 



where the two results correspond to the B s — > fo form 
factor evaluated at the leading order or not. Af'~^* 



are among the B s 



1 transition form factors and 



are taken from Rcf. 371. In Rcf. 



the ratio R^/^, can be inferred from the ratio of D s decay 



widths to /o7r + and </>7r + , obtaining Rf"/^ '- 
which is compatible with our result (|3"9"1) . 
Another relation has been also proposed in 



it was suggested that 
D s decay 

0.2 - 0.3, 



connect- 



ing Rfo/<p to a different observable in D s decays: 



R fo/<t> 



R fo/<l> 



(40) 

For this quantity the CLEO Collaboration has recently 
provided a measurement: Rf*/^ = (0.42 ± 0.11) [10| 
which is larger than our (|39j) . 

All the above considerations show that the mode B s — > 
J/t(>fo must be used, together with the golden mode 
B s — > J/ifxf), to measure the B s mixing phase, mainly 
because it provides us with a large number of events and 
does not require an angular analysis to separate different 
CP components of the final state. This is also the case of 
modes in which J / tp is replaced by a spin charmonium 
state, such as XcO, modulo the difficulty of the XcO re- 
construction. B s — > Xco4> will provide a side-check when 
the number of accumulated data will increase. Although 
B s — > Xco4> is a suppressed channel in naive factoriza- 
tion, its branching fraction may not be small due to the 
intermediate rescattering mechanism f38j or because of 



the contribution of nonfactorizablc diagrams [39j as in 
B — » XcoK. Analogously for B s — s- XcO(/>, the branching 
ratio of B — »• XcoK* has been measured [ioj . 



££(5° -> XcoA* ) 
KR(B- -> Xco A*") 



= (1.7 ±0.3 ±0.2) x 10" 4 , 
= (1.4 ±0.5 ±0.2) x 10~ 4 
< 2.1 x 10" 4 (90% CL) (41) 



and, on the basis of SU(?>)f symmetry, the branching 
fraction of B s — > Xco^ should be similar. 



0.6 r 

0.5 
0.4 

F!(0) 0.3 
0.2 
0.1 
0.0 L 




M 2 (GeV 2 ) 

Figure 3: Dependence of the D s — > fo form factors at q 2 = 
Fi(0) = F (0) on the Borel parameter M 2 . 
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V. DECAY D s ->■ f e + v 

By a suitable change of parameters in the sum rules in 
Section mi also the D s — s- f form factors can be com- 
puted and the branching ratio of the scmilcptonic decay 
D s — ^ /o e+jy can be predicted. We use m c — 1.4 GeV 
and t(D s ) = 0.5 ps [13; the threshold parameter is fixed 
to Sq s = (6.5 ± 1.0) GeV 2 . For the D s decay constant 
we use the value quoted by the Heavy Flavor Averaging 
Group: f Ds = (256.9 ± 6.8) MeV gH . The Borcl pa- 
rameter can be fixed requiring stability of the sum rule 
result with respect to M 2 variations. In Fig. |3] we plot 
versus M 2 ; the stability window is selected in 
the range M 2 = (5 ± 1) GeV 2 . We find 

F^ fo {0) = F^f°(0) = 0.30 ±0.03. (42) 

The q 2 dependence of the two form factors is displayed 
in Fig. |U The value of (|42|) is much smaller than in the 

1.0 p — . — . — . — i — . — . — . — i — . — . — . — i — . — . — . — i — . — . — . — , 

0.8 - 
0.6 - 

FM 1 ) \ 

0.4 - 



0.2 - 

o.o L — . — . — i — i — i — i — . — i — . — . — i — i — . — i — . — i — i — . — i — 

0.0 0.2 0.4 0.6 0.8 1.0 

q 2 (GeV 2 ) 

0.8 p — . — . — . — , — . — . — . — i — . — . — . — , — . — . — . — i — . — . — . — , 
0.6 - 




0.2 - 

0.0 Lj — ' — ' — ' — i — ' — ' — ' — i — ' — ' — ' — i — ' — ' — ' — i — ' — ' — ' — 
0.0 0.2 0.4 0.6 0.8 1.0 

9 2 (GeV 2 ) 

Figure 4: q 2 dependence of the D s — > fo form factors. 

D — > K case, for which the light-cone sum rule predic- 
tion is: Ff-> K (0) = 0.75±£og g|. We can understand 
this difference noticing that contribution of the fo twist-2 
LCDA in D s — > fo transition is small due to the differ- 
ent shape of the twist-2 fo distribution amplitude with 
respect to the case of K, The two LCDA are plotted in 
Fig.[5l where the position of the parameter uo, defined in 
Eq. (f2~Tj) , is also displayed (upper panel). The situation 
can be compared to the B s — > fo case, shown in the lower 
panel of the figure. Since the LCDA is integrated in the 
range [uo, 1], one can see that, in the case of Jo (980), the 
integral of the distribution amplitude gets two opposite 
contributions which tend to cancel each other, due to the 



presence of a zero in the DA. The zero is not present in 
the kaon DA, so that the integrated DA gives a much 
larger contribution. In the case of B s , the position of 
the parameter uq is such that the zero of the DA is not 
included in the integration region, so that no sizable dif- 
ference is expected between the fo and the kaon cases. 
This argument explains also why, compared with the re- 



1.5 




-1.5 u 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 u 

0.0 0.2 0.4 0.6 0.8 1.0 



U 



1.5 




-1.5 u — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — — 1 — 1 — 1 — 1 — 1 — u 
0.0 0.2 0.4 0.6 0.8 1.0 



U 

Figure 5: Shape of the twist-2 LCDA: -$ /o = -6u(l - 
u)B 1 Cl /2 (2u - 1) (dashed) and $ K (solid) taken from 
In the upper panel, the red line denotes the position of 
u® s — 0.334 fixed for the D s — > fo transition, while in 
the lower panel the red line corresponds to the position of 
Uq s — 0.684 at q 2 — in B s — > fo transition. 

suits of other approaches, our outcome are smaller. This 
can be noticed in Table IIV1 where we compare our re- 
sults for the D s — > fo form factors with other estimates 

The form factor F\(q 2 ) enters in the expression of the 



Table IV: D s -> /o(980) form factor at q 2 = 0, together 
with the results obtained by CLFD/DR [2{J, QCDSR 
and CLFQM [HI approaches. 



CLFD/DR 


QCDSR 


CLFQM This work 


Fi(0) 0.45/0.46 a 


1.7(0.27 ±0.02) 6 


0.434 0.30 ±0.03 



"using f Ds = 0.274 GeV 

'using fjj = 0.22±0.02 GeV; by using different input parameters 
two resuits are obtained, the first one in parentheses, the second 
one 1.7 times iarger. 
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differential decay rate 



dT( D s -> / e+„) G 2 F V 2 s \ 3 / 2 {m 2 Ds , m 2 Q ,q 2 ) 



2\|2 



dq 2 



192m 3 Ds ir 3 



(43) 



where the lepton mass is neglected. Since in D s — > foe + v 
the kincmatically accessible q 2 range is limited, the ap- 
plicable region for LCSR is narrow. One can fit the form 
factors in the spacelike region, for example — 2GeV 2 < 
q 2 < 0, and then extrapolate to the timelike region. How- 
ever, the result of the extrapolation strongly depends on 
the choice of the fitting region. Moreover, looking at Fig. 
HI one can notice that the q 2 dependence of F\ and Fq 
is mild. In view of this, we use a constant form fac- 
tor Fi(q 2 ) = Fi(0) to compute the branching ratio of 
D s /o e+l/ ; the result varies less than 10% including 
the q 2 dependence according to different fitting formulae. 
The obtained branching fraction is 



SK(£> a ->■ f e+v) = (2.0±g;t) x 10" 



(44) 



The modification due to radiative corrections can be es- 
timated as in the case of B s — » fo, finding F-j * - *' (0) = 
0.29 ' " 



+0.05 
-0.04- 



Let us consider the available experimental data. 
The CLEO Collaboration has measured the product of 
branching fractions [l(J 



/o(980)e+^) x £ft(/ -> 7T+7T-) = 
(0.20 ±0.03 ±0.01) x 10~ 2 , 



!B^(L» ; 

updating a previous determination Q 



(0.13 ±0.04 ±0.01) x 10~ 2 . 



(45) 



(46) 



Using experimental data provided by the BES Collabo- 
ration studying the processes \c0 — > fo (980) fo (980) — > 



iv+n-K+K- and X c0 /o(980)/ (980) 
K+K~K+K~ [H], CLEO quotes 



(47) 



S3?(/o -> tt+tt") = (50+£) x 10~ 2 
which, combined with (|45[) . gives 

■B^(L> S -> / e+i/) = (4.0 ± 0.6 ± 0.6) x 10~ 3 , (48) 
marginally consistent with our (|44l) . 



VI. CONCLUSIONS 

We have computed the B s — > fo transition form factors 
using light-cone QCD sum rules at leading order in the 
strong coupling constant, and also estimating the size 
of NLO corrections. The resulting form factors permit 
to predict the rates of B s — > fo£ + £~ and B s — » fovv 
decays, finding branching ratios accessible at future ma- 
chines, like a Super B factory, and at the LHCb exper- 
iment at CERN. The branching ratio of B s — > J/ipfo 
can be predicted under the factorization assumption: we 
find T>Ji(B s -> J/il)f )/'B3l(B a -> J/ip4>) - 0.2 - 0.3, 
thus the B s — > J/ipfo channel can be considered another 
promising mode to access the B s — B s mixing phase. We 
have also investigated the D s — > foe + v decay channel by 
the LCSR approach and compared the results to recent 
measurements . 
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